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The analysis of the radio light curves of the blazar CTA 102 during its 2006 flare revealed 
a possible interaction between a standing shock wave and a traveling one. In order to 
better understand this highly non-linear process, we used a relativistic hydrodynamic 
code to simulate the high energy interaction and its related emission. The calculated 
synchrotron emission from these simulations showed an increase in turnover flux density, 
Sm , and turnover frequency, Um , during the interaction and decrease to its initial values 
after the passage of the traveling shock wave. 
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1. Introduction 

The blazar CTA 102 (z=1.036) underwent a major flux density outburst in April 
2006 reaching a historical value of 10 Jy at 37 ghA The analysis of the single 
dish light curves revealed an unexpected double hump structure in the turnover 
frequency - turnover flux density [um — Sm) plane. A detailed modeling of this 
flare with the shock-in-jet model suggested a possible traveling-standing shock 
interaction as the origin of the flux density increas^. This scenario is supported 
by studies of the 15 GHz and 43 GHz parsec-scale morphological changes in the 
source from very-long-baseline interferometry (VLBI) images. These studies show 
a stationary feature at a distance of 0.2 mas (36 pc de-projected) away from the 
radio core (taken as the brightest feature at the base of the jet) during the period 
2003 -2011. At the time of the 2006 flare, feature ejected from the jet base passed 
through this stationary one, increased the observed flux densities and dragged the 
component outwards^. Standing shocks or re-coUimation shocks arise naturally in 
over-pressured jets and are characterized by a local increase in pressure, density and 
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magnetic fielciSl. Relativistic hydrodynamical simulations of superluminal sources 
have shown that the interaction of the standing shock with a travehng one could 
lead to the above mentioned observational effectiSl. 



2. Simulation 

For our simulation we used a relativistic hydrodynamic co based on a high- 
resolution shock-capturing scheme for relativistic hydrodynamicJ^. We simulated a 
relativistic jet consisting of relativistic electrons and sub-relativistic protons using 
an ideal gas equation of state with a constant adiabatic exponent 7 ~ 13/9. The 
initial jet radius was set to Rj = 0.3 pc, with velocity Vj = 0.99652 c, Mach number 
Mj = 3, density pj = 3.34 x 10~^^g/cm^ and pressure pj = 6.0 x 10~^ dyne/cm^. 
The ambient medium is assumed to be constant along the jet with density Pa = 
1.67 • 10~^'^g/cm'^ and pressure pj = 2.0 ■ 10^^ dyne/cm^, which leads to an over- 
pressure of the jet of dk = Pj/Pa = 3. The numerical resolution was set to 16 
cells per jet radii. This set of parameters corresponds to a jet with a kinetic power 
Lj - 3.0 X lO''^ erg/s, which is a typical value for FR II galaxies. After the jet 
reached a steady state, we injected a perturbation with pressure and density values 
four times larger than the values of the jet. For the duration of the injection event we 
assumed 3 months which was comparable to the observed duration of the increase 
of the flux density in CTA 102 during the 2006 radio flare. From the pressure, 
density, and fluid velocity distributions obtained, we calculated the synchrotron 
emission following refs. |8l [9] and references therein, assuming that the magnetic 
energy density is a fraction =0.1 of the equipartition magnetic field and that a 
fraction eg = 0.1 of the shock energy is used to accelerate the non-thermal particles. 
The results of this simulation are presented in Fig. [1] and will be described in the 
following section. 



3. Results 

Due to the initial over-pressure in the jet, the fluid expanded radially after leaving 
the jet nozzle. During this process, a rarefied region forms close to the jet bound- 
aries. In the inner jet regions, however, a conical shock from the injection generates 
a Mach disk around the jet axis at a distance of z 10 pc from the injection. The 
final structure could be best described by an inner spine (r < 0.3 pc), characterized 
by density, pressure, and velocity values comparable to the ones at the jet nozzle, 
and an outer sheath (0.3 pc < r < 1 pc) with decreasing density and pressure profiles 
with distance and increasing fluid velocity, as expected from an adiabatic expan- 
sion. Downstream of the Mach disk, the fluid expands conically, showing a shallow 
acceleration and creating the first re-collimation shock at a distance of z ~ 33 pc 
from the jet nozzle. At this point, the fluid decelerates and the density and pressure 
values jump nearly to their initial values (factor 15 in pressure and factor 20 in den- 
sity). The outer sheath shows a similar behaviour; the radial expansion reaches a 
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maximum around a distance of z 25 pc and starts to recollimate. This collimation 
gives rise to inward traveling sound waves which interact with the waves created by 
the colhmation of the spine and form an extended region of increased pressure and 
density between 33pc < z < 50 pc (see top panel in Fig. [Ij. As mentioned before, 
we injected a perturbation after the jet reached a steady state. In the following, we 
concentrate on the variations of pressure, density and velocity of the inner spine. 
The injected perturbation travels along a decreasing pressure and density profile 
producing an increase by a factor 4 in pressure and density (the jumps in pressure 
and density were decreasing with distance from the nozzle). However, at the posi- 
tion of the first re-collimation shock {z ^ 33 pc), there is a partial reflection of the 
perturbation and it can not pass through the standing shock at once (middle panel 
in Fig. [1]). Instead, we found that the perturbation piles up to a certain pressure 
and density value before it can cross the re-collimation shock. The piled up material 
leads to an increase of a factor 12 in pressure and a factor of 6 in density at the 
location of the re-collimation shock, with respect to the steady-state values (bottom 
panel in Fig.[l]). At the same time, the standing shock is dragged downstream for a 
specific time before it is re-established at its initial position. We calculated the ob- 
served non-thermal emission using the known distance to the source, a viewing angle 
of 2.6° and a spectral slope, s = 2.5. For the steady state jet we obtain a turnover 
frequency of Vm ^ 1 GHz and a turnover flux density of 5*^ ^ 0.03 Jy at the position 
of the re-collimation shock and values of ~ 15 GHz and Sm ^ 0.24 Jy during the 
interaction with the traveling shock. After the passage of the traveling shock wave 
through the standing one, the emission values decreased back to the ones obtained 
for the steady state. 



4. Conclusions 

Our simulation of the interaction between a traveling shock and a stationary shock 
showed that this process can lead to a strong increase in the observed non-thermal 
radiation. However, the turnover values obtained during the shock-shock interaction 
are lower than the observed ones in CTA 102 during the 2006 outburst (i^™ ~ 
50 GHz, Sm ~ 8 Jy). This could be an indication for: 

(1) The density in relativistic jets is larger than 20% that of the ambient medium 
(or the ambient medium is denser). 

(2) The influence of the magnetic field is higher than assumed i. e., ef, > 0.1. 

(3) Shocks are more efficient in accelerating particles than assumed, i. e., eg > 0.1. 

To verify the assumptions above, two directions of research are planned. We will per- 
form a deep analysis of the existing multi-frequency VLBI observations of CTA 102 
and an extended parameter study for the simulated shock-shock interactions. 
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Fig. 1. Results of our RHD simulation of the interaction between a traveling and a stationary 
shock in a rclativistic jet. Shown are the logarithm of the rest mass density for the steady-state jet 
(top panel), partial reflection of the traveling shock wave at the onset of the re-coUimation shock 
(middle panel) and the crossing of re-coUimation shock by the traveling shock (bottom panel) 
Notice the increased density value during the shock-shock interaction. 
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